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Abstract. We show that the spin-and-space unrestricted Hartree-Fock method, in conjunction with the 
companion step of the restoration of spin and space symmetries via Projection Techniques (when such 
symmetries are broken), is able to describe the full range of couplings in two-dimensional double quantum 
dots, from the strong-coupling regime exhibiting delocalized molecular orbitals to the weak-coupling and 
dissociation regimes associated with a Generalized Valence Bond combination of atomic-type orbitals 
localized on the individual dots. The weak-coupling regime is always accompanied by an antiferromagnetic 
ordering of the spins of the individual dots. The cases of dihydrogen (Hh, 2e) and dilithium (Li2, 6e) 
quantum dot molecules are discussed in detail. 

PACS. 73. 21. La Quantum dots - 85.35.-p Nanoelectronic devices - 31.15.Rh Valence bond calculations 



1 Introduction 



Two-dimensional (2D) Quantum Dots (QD's) are usually 
referred to as artificial atoms, a term suggestive of strong 
similarities between these manmade devices and the phys- 
ical behavior of natural atoms. As a result, in the last few 
years, an intensive theoretical effort@,|||,|,|,|,0|,|,|lO| 
has been devoted towards the elucidation of the appropri- 
ate analogies and/or differences. Recently, we have shown 
|0,||g] that, even in the absence of a magnetic field, the 
most promising analogies are found mainly outside the 
confines of the central-field approximation underlying the 
Independent-Particle Model (IPM) and the ensuing phys- 
ical picture of electronic shells and the Aufbau Princi- 
ple. Indeed, as a result of the lower electronic densities 
in QD's, strong e — e correlations can lead (as a func- 
tion of the ratio Rw between the interelectron repulsion 
and the zero-point kinetic energy) to a drastically differ- 
ent physical regime, where the electrons become localized, 
arranging themselves in concentric geometric shells and 
forming electron molecules. In this context, it was found 
[g,U that the proper analogy for the particular case of 
a 2e QD is the collective-motion picture reminiscent of 
the fleeting and rather exotic phenomena of the doubly- 
excited natural helium atom, where the emergence of a 
"floppy" trimeric molecule (consisting of the two localized 
electrons and the heavy a-particle nucleus) has been well 
established @|2). 

A natural extension of this theoretical effort has also 
developed in the direction of lateral 2D QD Molecules 
(QDM's, often referred to as artificial molecules), aiming 
at elucidating M,nfl,0,n4,15,lq] the analogies and differ- 
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ences between such artificially fabricated nanostrustures 
and the natural molecules. 

In this paper, we address the interplay of coupling and 
dissociation in QDM's, and its effects concerning the ap- 
pearance of ferromagnetic versus antiferromagnetic order- 
ing. We will show that this interplay relates directly to the 
nature of the coupling in the artificial molecules, and in 
particular to the question whether such coupling can be 
descibed by the Molecular Orbital (MO) Theory or the 
Valence Bond (VB) Theory in analogy with the chemical 
bond in natural molecules. 

Furthermore, we show that the onset at a moderate 
interdot barrier or interdot distance do, as well as the per- 
manency for all separations d > do, of spontaneous magne- 
tization and ferromagnetic ordering predicted for Double 
QD's by local-spin-density (LSD) calculations [^3||lJ] is 
an artifact of the MO structure implicit in the framework 
of these density-functional calculations. 

We utilize a self-consistent-field theory which can go 
beyond the MO approximation, namely the spin-and-space 
unrestricted Hartree-Fock (sS-UHF) , which was introduced 
by us M,p| for the description of the many-body prob- 
lem of both single f||| and molecular @ QD's. This sS- 
UHF employs N (where N is the number of electrons) 
orbital-dependent, effective (mean- field) potentials and it 
differs from the more familiar [[l7j restricted HF (RHF) in 
two ways: (i) it relaxes the double- occupancy requirement, 
namely, it employs different spatial orbitals for the two dif- 
ferent (i.e., the up and down) spin directions [DODS, thus 
the designation "spin (s) unresricted" ] , and (ii) it relaxes 
the requirement that the electron orbitals be constrained 
by the symmetry of the external confining field [thus the 
designation "space (S) unrestricted"]. 
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In this paper we will show that the solutions with bro- 
ken space symmetry allowed in QDM's by the sS-UHF 
provide a natural vehicle for formulating a Generalized 
Valence Bond (GVB) theory. Furthermore, they result in 
an antiferromagnetic ordering of the molecular ground 
state, in contrast to the ferromagnetic ordering of the 
MO method which is associated with HF solutions that 
preserve the space symmetry, namely those derived from 
the fully restricted Hartree-Fock (RHF) or the spin un- 
restricted (but not space unrestricted) Hartree-Fock (s- 
UHF) approaches. 



(v+0.5)hu> y i + hi denotes the y- eigenvalues. The matching 
conditions at y = for the left and right domains yield the 
^-eigenvalues and the eigenfunctions Y v (y) (m is integer 
and v is in general real) . 

In this paper, we will focus on the zero-field case (B = 
0) and we will limit ourselves to symmetric (homopolar) 
QDM's, i.e., hu> x = hu> y i = hu> y 2 = hojQ, with equal well- 
depths of the left and right dots, i.e., hi = h 2 = 0. In all 
cases, we will use Tiujq = 5 meV and m* = 0.067m e (this 
effective-mass value corresponds to GaAs). 



2 The two-center-oscillator confining 
potential 

In the 2D two-center-oscillator (TCO), the single-particle 
levels associated with the confining potential of the artifi- 
cial molecule are determined by the single-particle hamil- 
tonian [frl 
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y k = y — y k with k = 1 for y < (left) and k = 2 
for y > (right), and the h k s control the relative well- 
depth, thus allowing studies of hetero-QDM's. x denotes 
the coordinate perpendicular to the interdot axis (y). T — 
(p-eA/c) 2 /2m*, with A = 0.5(- By, Bx,0), and the last 
term in Eq. (HI) is the Zeeman interaction with g* being 
the effective g factor and fis the Bohr magneton. Here we 
limit ourselves to systems with hui x i = hw x i = h~to x - The 
most general shapes described by H are two semiellipses 
connected by a smooth neck [V ne ck(y)}- J/i < and y 2 > 
are the centers of these semiellipses, d = y 2 — Vi is the 
interdot distance, and m* is the effective electron mass. 

For the smooth neck, we use V neck {y) = \m* '^l k [c k y k '+ 
dky' k 4 ]0{\y\ - \y k \), where 6{u) = for u > and 9{u) = I 
for u < 0. The four constants c k and dk can be expressed 



via two parameters, as follows: (- 
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and 



d k = (1 — 3e b k )/yl, where the barrier-control parameters 
e fe = (Vb — h k )/Vo k are related to the actual (controlable) 
height of the bare barrier (Vf,) between the two QD's, and 
Vok = m*wl k yl/2 (for h x = h 2 , V 01 =V 02 =V Q ). 

The single-particle levels of H, including an external 
perpendicular magnetic field B, are obtained by numer- 
ical diagonalization in a (variable-with-separation) basis 
consisting of the eigenstates of the auxiliary hamiltonian: 



H = 
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This eigenvalue problem is separable in x and y, i.e., the 
wave functions are written as <t> mu (x,y) — X m (x)Y v (y). 
The solutions for X m (x) are those of a one-dimensional 
oscillator, and for Y v (y) they can be expressed through 
the parabolic cylinder functions Eq] U[a k , ( — l) fe Cfc] , where 
^ fe = y' k ^/2m*L0 yk /h, a k = {—E y + h k )/(huo yk ), and E y = 



3 The Many-Body Hamiltonian 

The many-body hamiltonian 7i for a dimcric QDM com- 
prising N electrons can be expressed as a sum of the single- 
particle part H{i) defined in Eq. (fil) and the two-particle 
interelectron Coulomb repulsion, 
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where k is the dielectric constant and r^ denotes the rel- 
ative distance between the i and j electrons. 

As we mentioned in the introduction, we will use the 
sS-UHF method for determining at a first level an approx- 
imate solution of the many-body problem specified by the 
hamiltonian <M). The sS-UHF equations were solved in the 
Pople-Nesbet-Roothaan formalism pj| using the interdot- 
distance adjustable basis formed with the eigenfunctions 
^mv{x,y) of the TCO defined in section 2. 

As we will explicitly illustrate in section 5 for the case 
of the H2-QDM, the next step in improving the sS-UHF 
solution involves the use of Projection Techniques in rela- 
tion to the UHF single Slater determinant. 



4 A first example of a homopolar 
two-dimensional artificial molecule: the 
lateral Li 2 -QDM 

As an illustrative example, we consider here an open-shell 
(N = 6 electrons) QDM made of two QD's (hence the 
name L12-QDM), with an interdot distance of d = 70 nm 
and an interdot barrier height of V& = 10 meV. The value 
of the dielectric constant is first taken to be k = 20.0. This 
value is higher than the value «(GaAs)= 12.9 for GaAs 
and may be viewed as resulting from screening produced 
by the external charges residing in the gates and/or the 
finite height (in the z direction) of the dot. This value of n 
is chosen such that the e — e repulsion is not strong enough 
to precipitate individual electron localization [p0| within 
each individual dot (that is, to inhibit Wigner crystalliza- 
tion on each of the dots, see Ref. M). This example thus 
conforms better to the case of a natural Li2 molecule. It 
was earlier presented briefly in Ref. H in connection with 
the formation of electron puddles (EP's). Here we will 
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FIGURE 1 



Fig. 1. Lateral Li2-QDM: s-UHF spin-up occupied orbitals 
(modulus square) and total charge (CD) and spin (SD) densi- 
ties (top row) for the P — 2 spin polarized case. The numbers 
displayed with each orbital are their s-UHF eigenenergies in 
meV. The two spin-down orbitals are not displayed; they are 
similar to the o" s (T) and the a u (T) an d have energies of 17.10 
meV and 17.27 meV, respectively. The number displayed with 
the total charge density (top leftO is the s-UHF total energy in 
meV. Distances are in nm and the electron densities in 1CF 4 



nm . The choice of parameters is: m* 
meV, d = 70 nm, V b = 10 meV, k = 20. 



0.067m e , huio 



FIGURE 2 



Fig. 2. Lateral Li2 QDM: RHF occupied orbitals (modulus 
square) and charge (CD) and spin (SD) densities (top row) for 
the P — spin unpolarized case. The numbers displayed with 
each orbital are their RHF eigenenergies in meV, while the up 
and down arrows indicate electrons with up or down spin, re- 
spectively. The number displayed with the total charge density 
is the RHF total energy in meV. Distances are in nm and the 
electron densities in 10~ 4 nm" 2 . The choice of parameters is 
the same as in Fig. 1. 



present a detailed study of it by contrasting the descrip- 
tions resulting from both the space-symmetry preserving 
(RHF or s-UHF) and non-preserving (sS-UHF) methods. 
We start by presenting HF results which preserve the 
space symmetry of the QDM; the results corresponding to 



the P = N 7; — N 1= 2 and the P = spin polarizations 
are displayed in Fig. 1 and Fig. 2, respectively. The P = 2 
case is an open-shell case and is treated within the (spin 
unrestricted) s-UHF method, in analogy with the standard 
practice for open-shell configurations in Quantum Chem- 
istry pH ; the P = case is of a closed-shell- type and is 
treated within the RHF approach. 

The corresponding energies are E s uhf (P = 2) = 74.137 
meV and Erhf{P = 0) = 75.515 meV, i.e., the space- 
symmetry preserving HF variants predict that the polar- 
ized state is the ground state of the molecule. Indeed, 
by preserving the symmetries of the confining potential, 
the associated orbitals are clearly of the MO-type. In this 
MO picture the L12-QDM exhibits an open shell struc- 
ture (for 2D QDM's, the closed shells correspond to N = 
4, 12, 24, ... electrons, that is shell closures at twice the val- 
ues corresponding to an individual harmonic 2D QD) and 
thus Hund's first rule should apply for the two electrons 
outside the core closed shell of the first 4 electrons. 

Fig. 1 displays for the P — 2 spin polarization the 
four occupied spin-up s-UHF orbitals (a g , a u , ir x ,gj 7tx,«); 
as well as (top of Fig. 1) the total charge density (CD, 
sum of the N | +N I electron densities) and spin density 
(SD, difference of the N f — N J. electron densities). The 
two occupied spin-down orbitals (cr' g , <j' u ) are not shown: 
in conforming with the DODS approximation, they are 
slightly different from the a g {\) and <j u (1) ones. 

Fig. 2 displays the corresponding quantities for the 
P = case. In both (the P = 2 and P = 0) cases, the 
orbitals are clearly molecular orbitals delocalized over the 
whole QDM. The P = 2 polarization (Fig. 1) exhibits the 
molecular configuration <Jg<j' (T„<r„ 7t^ g 7r^ u ; this configu- 
ration changes to u^a^K^ in the P — case (Fig. 2). 

Fig. 3 displays the corresponding quantities for the 
P = state calculated using the (spin-and-space unre- 
stricted) sS-UHF. This state exhibits a breaking of space 
symmetry (the reflection symmetry between the left and 
right dot). Unlike the MO's of the P = RHF solution 
(Fig. 2), the sS-UHF orbitals in Fig. 3 are well localized 
on either one of the two individual QD's and strongly re- 
semble the atomic orbitals (AO's) of an individual Li-QD, 
i.e., they are of Is and lp x type. Comparing with the MO 
case of Fig. 2, one sees that the symmetry breaking did 
not greatly influence the total charge density. However, 
a dramatic change appeared regarding the spin densities. 
Indeed the SD of the sS-UHF solution (see top right of 
Fig. 3) exhibits a prominent spin density wave (SDW) as- 
sociated with an antiferromagnetic ordering of the coupled 
individual QD's. Formation of such SDW ground states in 
QDM's is accompanied by electron (orbital) localization 
on the individual dots, and thus in Ref. |j we proposed 
for them the name of Electron Puddles (EP's). Notice that 
the EP's represent a separate class of symmetry broken 
solutions, different M from the SDW's that can develop 
within a single QD p2| and whose formation involves the 
relative rotation of delocalized open-shell orbitals, instead 
of electron localization. 

The sS-UHF total energy for the P — unpolarized 
case is E s suhf{P = 0) = 74.136 meV, i.e., the sym- 
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FIGURE 3 



FIGURE 4 



Fig. 3. Lateral Li2-QDM: sS-UHF occupied orbitals (modulus 
square) and charge (CD) and spin (SD) densities (top row) for 
the P = spin unpolarized case. The numbers displayed with 
each orbital are their sS-UHF eigenenergies in meV, while the 
up and down arrows indicate an electron with an up or down 
spin. The number displayed with the charge density is the sS- 
UHF total energy in meV. Distances are in nm and the electron 



densities in 10 
as in Fig. 1. 



nm . The choice of parameters is the same 



mctry breaking produces a remarkable gain in energy of 
1.379 meV. As a result, the unpolarized state is the ground 
state, while the ferromagnetic ordering predicted by the 
RHF is revealed to be simply an artifact of the MO struc- 
ture implicit in this level of approximation. Notice that 
the symmetry-broken unpolarized state is only 0.001 meV 
lower in energy than the P — 2 polarized one, namely the 
two states are practically degenerate, which implies that 
for the set of parameters employed here the QD molecule 
is located well in the dissociation regime. 

Although the symmetry breaking within the HF theory 
allows one to correct for the artifact of the spontaneous 
polarization exhibited by the MO approaches, i.e., the the 
fact that the MO polarized solutions are the lowest in en- 
ergy upon separation, the resulting sS-UHF many-body 
wave function violates both the total spin and the space 
reflection symmetries. Intuitively further progress can be 
achieved as follows: the orbitals (see Fig. 3) of the sS-UHF 
can be viewed as optimized atomic orbitals (OAO's) to 
be used in constructing a many-body wave function con- 
forming to a Generalized Valence Bond (GVB) structure 
and exhibiting the correct symmetry properties. Starting 
with the sS-UHF symmetry-breaking Slater determinant, 
this program can be systematically carried out within the 
framework of the theory of Restoration of Symmetry (RS) 
via the so-called Projection Techniques. Instead of the 
more complicated case of the Li2-QDM, and for reasons 
of simplicity and conceptual clarity, we will present below 
the case of the QD molecular hydrogen as an illustrative 
example of this RS procedure. 



Fig. 4. Lateral H2 QDM: Occupied orbitals (modulus square, 
bottom half) and total charge (CD) and spin (SD) densities 
(top half) for the P — spin unpolarized case. Left column: 
RHF results. Right column: sS-UHF results exhibiting a break- 
ing of the space symmetry. The numbers displayed with each 
orbital are their eigenenergies in meV, while the up and down 
arrows indicate an electron with an up or down spin. The num- 
bers displayed with the charge densities are the total energies 
in meV. Unlike the RHF case, the spin density of the sS-UHF 
exhibits a well developed spin density wave. Distances are in 



parameters is: m* = 0.067m e 
V b = 4.95 meV, k = 20. 



Tiuio = 5 meV, d = 30 nm, 



5 Artificial molecular hydrogen (H 2 -QDM) in 
a Generalized Valence Bond Approach 

5.1 The sS-UHF description 

We turn now to the case of the H2-QDM. Fig. 4 displays 
the RHF and sS-UHF results for the P = case (sin- 
glet) and for an interdot distance d = 30 nm and a bar- 
rier V b = 4.95 meV. In the RHF (Fig. 4, left), both the 
spin-up and spin-down electrons occupy the same bonding 
(<j g ) molecular orbital. In contrast, in the sS-UHF result 
the spin-up electron occupies an optimized Is (or Isl) 
atomic- like orbital (AO) in the left QD, while the spin 
down electron occupies the corresponding Is' (or Isr) 
AO in the right QD. Concerning the total energies, the 
RHF yields E RHF (P = 0) = 13.68 meV, while the sS- 
UHF energy is E s suhf(P = 0) = 12.83 representing a 
gain in energy of 0.85 meV. Since the energy of the triplet 
is E SUHF {P'= 1) = E sSUHF (P = 1) = 13.01 meV, the 
sS-UHF singlet conforms to the requirement that for two 
electrons at zero magnetic field the singlet is always the 
ground state; on the other hand the RHF MO solution 
fails in this respect. 



5.2 Projected wave function and restoration of the 
broken symmetry 

To make further progress, we utilize the spin projection 
technique to restore the broken symmetry of the sS-UHF 
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determinant (henceforth we will drop the prefix sS when 
referring to the sS-UHF determinant), 



V2&uhf(1,2) = 



«(r x )o(l) v( ri )0(l) 
«(r 2 )a(2) v(r 2 )/3(2) 

= \u{l)v(2)> , 



(4) 



where w(r) and v(r) are the Is (left) and Is' (right) lo- 
calized orbitals of the sS-UHF solution displayed in the 
right column of Fig. 4; a and [3 denote the up and down 
spins, respectively. In Eq. (0) we also define a compact 
notation for the ^uhf determinant, where a bar over a 
space orbital denotes a spin-down electron; absence of a 
bar denotes a spin-up electron. 

&uhf(1, 2) is an eigenstate of the projection S z of the 
total spin S = si + s 2 , but not of S 2 . One can generate 
a many-body wave function which is an eigenstate of S 2 
with eigenvalue s(s + 1) by applying the following projec- 
tion operator introduced by Lowdin [f23|,p4|, 



Ps^U 



s't^s 



S 2 - s'{s' + l)h 2 
Ks + l)-s'(s' + l)]^ 2 



(5) 



where the index s' runs over the eigenvalues of S 2 . 

The result of S 2 on any UHF determinant can be cal- 
culated with the help of the expression, 



s 2 <Puhf = n 2 



{n a -np) 2 /A + n/2 + 






*. 



UHF 



(6) 

where the operator wij interchanges the spins of electrons 
i and j provided that their spins are different; n a and np 
denote the number of spin-up and spin-down electrons, 
respectively, while n denotes the total number of electrons. 
For the singlet state of two electrons, n a = np = 1, 
n = 2, and S 2 has only the two eigenvalues s = and 
s = 1. As a result, 

2y/2Pb$ UHF (l,2) = (l-C7i2)>/MW(l,2) 

= ]«(l)v(2) > -|fi(l)v(2) > . (7) 

In contrast to the single-determinantal wave functions of 
the RHF and sS-UHF methods, the projected many-body 
wave function (0) is a linear superposition of two Slater de- 
terminants, and thus represents a corrective step beyond 
the mean-field approximation. 

Expanding the determinants in Eq. (0), one finds the 
equivalent expression 

2P & UHF (1, 2) = (u(ri)»(r 2 ) + u(r 2 )v(ri)) X (0, 0) , (8) 

where the spin eigenfunction for the singlet is given by 

X (s = 0, S z = 0) = (a(l)/3(2) - a (2)/3(l))/V2 . (9) 

Eq. (0) has the form of a Heitler-London (HL) J25[ or 
valence bond |2(],|27],g8| wave function. However, unlike 
the HL scheme which uses the orbitals 4>l(?) and (j>R{r) 



of the separated (left and right) atoms [E9| , expression M) 
employs the sS-UHF orbitals which are self-consistently 
optimized for any separation d and potential barrier height 
Vb- As a result, expression (g) can be characterized as 
a Generalized Valence Bond (GVB) p9 wave function. 
Taking into account the normalization of the spatial part, 
we arrive at the following improved wave function for the 
singlet state exhibiting all the symmetries of the original 
many-body hamiltonian, 



« 



GVB 



T, 2) = N + V2P & UH f{1, 2) , (10) 



where the normalization constant is given by 



N + = l/^l + S 2 lv , (11) 

S uv being the overlap integral of the u(r) and v(r) orbitals, 

(12) 



(13) 



Suv = / u(r)v(r)dr 



The total energy of the GVB state is given by 



Egvb — N + [h uu + h t 



Zdl.rllfl 



UV""UV 



■J 11 



■K n 



where h is the single-particle part of the total hamiltonian 
(H), and J and K are the direct and exchange matrix el- 
ements associated with the e — e repulsion e 2 / KT\2- For 
comparison, we give also here the corresponding expres- 
sion for the HF total energy either in the RHF (v = u) or 
sS-UHF case, 



E 



HF 



<Ju 



(14) 



For the triplet, the projected wave function coincides 
with the original HF determinant, so that the correspond- 
ing energies in all three approximation levels are equal, 
i.e. 



E GVB 



e rhf 



E UHF- 



5.3 Comparison of RHF, sS-UHF and GVB results 

A major test for the suitability of different methods to the 
case of QDM's is their ability to properly describe the dis- 
sociation limit of the molecule. The H2-QDM dissociates 
into two non-interacting QD hydrogen atoms with arbi- 
trary spin orientation. As a result, the energy difference, 
As = E s — E t , between the singlet and the triplet states 
of the molecule must approach the zero value from below 
as the molecule dissociates. To theoretically generate such 
a dissociation process, we keep the separation d constant 
and vary the height of the interdot barrier Vf>; an increase 
in the value of VJ, reduces the coupling between the indi- 
vidual dots and for sufficiently high values we can always 
reach complete dissociation. 

Fig. 5 displays the evolution in zero magnetic field of 
As as a function of Vb and for all three approximation 
levels, i.e., the RHF (MO Theory, top line), the sS-UHF 
(middle line), and the GVB (bottom line). The interdot 
distance is the same as in Fig. 4, i.e., d = 30; the case 
of k — 20 is shown at the top panel, while the case of 
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Fig. 5. Lateral H2-QDM: The energy difference between the 
singlet and triplet states according to the RHF (MO Theory, 
top line), the sS-UHF (middle line), and the GVB approach 
(Projection Method, bottom line) as a function of the interdot 
barrier VJ,. For VJ, = 25 meV complete dissociation has been 
clearly reached. Top frame: k — 20. Bottom frame: k — 40. The 
choice of the remaining parameters is: m* = 0.067m,;, huo — 5 
meV and d = 30 nm. 



a weaker e — e repulsion is displayed for k — 40 at the 
bottom panel. 

We observe first that the MO approach fails completely 
to describe the dissociation of the HVQDM, since it pre- 
dicts a strongly stabilized ferromagnetic ordering in con- 
tradiction to the expected singlet-triplet degeneracy upon 
full separation of the individual dots. A second observation 
is that both the sS-UHF and the GVB solutions describe 
the dissociation limit (Ae — > for Vj, — > 00) rather well. 
In particular in both the sS-UHF and the GVB methods 
the singlet state remains the ground state for all values of 
the interdot barrier. Between the two singlets, the GVB 
one is always the lowest, and as a result, the GVB method 
presents an improvement over the sS-UHF method both 
at the level of symmetry preservation and the level of en- 
ergetics. 

It is interesting to note that for k = 40 (weaker e — e 
repulsion) the sS-UHF and GVB solutions collapse to the 
MO solution for smaller interdot-barrier values Vj, < 2.8 
meV. However, for the stronger e — e repulsion (k = 20) 
the sS-UHF and GVB solutions remain energetically well 
below the MO solution even for Vb = 0. Since the sepa- 



ration considered here (d = 30 nm) is a rather moderate 
one (compared to the value ^0 == 28.50 nm for the extent 
of the Is AO), we conclude that there is a large range of 
materials parameters and interdot distances for which the 
QDM's are weakly coupled and cannot be described by 
the MO theory. 



6 Conclusions 

We have shown that the sS-UHF method, in conjunction 
with the companion step of the restoration of symme- 
tries when such symmetries arc broken, is able to describe 
the full range of couplings in a QDM, from the strong- 
coupling regime exhibiting delocalized molecular orbitals 
to the weak-coupling one associated with Hcitler-London- 
type combinations of atomic orbitals. 

The breaking of space symmetry within the sS-UHF 
method is necessary in order to properly describe the weak- 
coupling and dissociation regimes of QDM's. The break- 
ing of the space symmetry produces optimized atomic-like 
orbitals localized on each individual dot. Further improve- 
ment is achieved with the help of Projection Techniques 
which restore the broken symmetries and yield multide- 
tcrminantal many-body wave functions. The method of 
the restoration of symmetry was explicitly illustrated for 
the case of the molecular hydrogen QD (H2-QDM, see sec- 
tion 5). It led to the introduction of a Generalized Valence 
Bond many-body wave function as the appropriate vehicle 
for the description of the weak-coupling and dissociation 
regimes of artificial molecules. In all instances the weak- 
coupling regime is accompanied by an antiferromagnetic 
ordering of the spins of the individual dots. 

Additionally, we showed that the RHF, whose orbitals 
preserve the space symmetries and are delocalized over the 
whole molecule, is naturally associated with the molecular 
orbital theory. As a result, and in analogy [pj||27|,|3l| with 
the natural molecules, it was found that the RHF fails 
to describe the weak-coupling and dissociation regimes 
of QDM's. It can further be concluded that the sponta- 
neous polarization and ferromagnetic ordering predicted 
for weakly-coupled double QD's in Rcfs. 0,0] is an arti- 
fact of the molecular-orbital structure implicit [B2| in the 
framework of LSD calculations. 
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The early electronic model of valence was primarily devel- 
oped by G.N. Lewis who introduced a symbolism where 
an electron was represented by a dot, e.g., H:H, with a 
dot between the atomic symbols denoting a shared elec- 
tron. Later in 1927 Heitler and London formulated the 
first quantum mechanical theory of the pair-electron bond 
for the case of the hydrogen molecule. The theory was sub- 
sequently developed by Pauling and others in the 1930's 
into the modern theory of the chemical bond called the 
Valence Bond Theory. 

Refs. [p_5j,[16j have studied, as a function of the magnetic 
field, the behavior of the singlet-triplet splitting of the QD 
H2 by diagonalizing the two-electron hamiltonian inside 
the minimal four-dimensional basis formed by the products 
0i(ri)0£(r 2 ), <J>l(ti)<J>r(t2), 0it(ri)0z,(r a ), </>fl(ri)0H(r 2 ) 
of the Is orbitals of the separated QD's. This Hubbard- 
type method Jl5| (as well as the refinement employed by 
Ref. [1m of enlarging the minimal two-electron basis to 
include the p orbitals of the separated QD's) is an im- 
provement over the simple HL method (see Ref. fl5|), but 
apparently it is only appropriate for the weak-coupling 
regime at sufficiently large distances and/or interdot bar- 
riers. Our method is free of such limitations, since we em- 
ploy an interdot-distance adjustable basis (see section 2) 
of at least 66 spatial TCO molecular orbitals when solv- 
ing for the sS-UHF ones. Even with consideration of the 
symmetries, this amounts to calculating a large number of 
two-body Coulomb matrix elements, of the order of 10 6 . 
More precisely our GVB method belongs to a class of Pro- 
jection Techniques known as Variation before Projection, 
unlike the familiar in chemistry GVB method of Goddard 
and coworkers [W.A. Goddard III et al., Ace. Chem. Res. 
6, 368 (1973)], which is a Variation after Projection [see 
P. Ring and P. Schuck, The Nuclear Many-Body Problem 
(Springer, New York, 1980) ch. 11]. 
see section 3.8.7 in Ref. Jl9[ . 

Symmetry breaking in coupled QD's within the LSD has 
been explored by J. Kolehmainen et al. [Eur. Phys. J. D 13, 
731 (2000)]. However, unlike the HF case for which a fully 
developed theory for the restoration of symmetries has long 
been established (see, e.g., the book by Ring and Schuck 
in Ref. [pOl), the breaking of space symmetry within the 
spin-dependent density functional theory poses a serious 
dilemma [J.P Perdew et al., Phys. Rev. A 51, 4531 (1995)]. 
This dilemma has not been fully resolved todate; several 
remedies (like Projection, ensembles, etc.) are being pro- 
posed, but none of them appears to be completely devoid 
of inconsistencies [A. Savin in Recent developments and ap- 
plications of modern density functional theory, edited by 
J.M. Seminario (Elsevier, Amsterdam, 1996) p. 327]. In 
addition, due to the unphysical self-interaction error, the 
density-functional theory is more resistant against sym- 
metry breaking [see R. Bauernschmitt and R. Ahlrichs, J. 
Chem. Phys. 104, 9047 (1996)] than the sS-UHF, and thus 
it fails to describe a whole class of broken symmetries in- 
volving electron localization, e.g., the formation at B — 
of Wigner molecules in QD's (see footnote 7 in Ref. JJ), 
the hole trapping at Al impurities in silica [J. Laegsgaard 
and K. Stokbro, Phys. Rev. Lett. 86, 2834 (2001)], or the 
interaction driven localization-delocalization transition in 
d- and /- electron systems, like Plutonium [SY. Savrasov 
et al, Nature 410, 793 (2001)]. 
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